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Abstract —We illustrate potential benefits of using massive 
antenna arrays for wireless energy transfer (WET). Specifically, 
we analyze probability of outage in WET over fading channels 
when a base station (BS) with multiple antennas beamforms 
energy to a wireless sensor node (WSN). Our analytical results 
show that by using massive antenna arrays, the range of WET 
can be increased for a given target outage probability. We prove 
that by using multiple-antenna arrays at the BS, a lower downlink 
energy is required to get the same outage performance, resulting 
into savings of radiated energy. We show that for energy levels 
used in WET, the outage performance with least-squares or 
minimum mean-square error channel estimates is same as that 
obtained based on perfect channel estimates. We observe that a 
strong line-of-sight component between the BS and WSN lowers 
outage probability. Furthermore, by deploying more antennas at 
the BS, a larger energy can be transferred reliably to the WSN 
at a given target outage performance for the sensor to be able 
to perform its main tasks. In our numerical examples, the RF 
power received at the input of the sensor is assumed to be on the 
order of a mW, such that the rectenna operates at an efficiency 
in the order of 50%. 

Index Terms —Wireless energy transfer, massive MIMO, beam¬ 
forming, outage probability, array gain 

I. Introduction 

Wireless energy transfer (WET) is a promising energy har¬ 
vesting technology where the destination node harvests energy 
from electromagnetic radiations instead of traditional wired 
energy sources 0. The use of WET can help increase the 
battery-lifetime of energy-constrained wireless sensor nodes 
(WSNs) that are used for applications such as intelligent 
transportation, intrusion detection, and aircraft structural moni¬ 
toring 0 . Furthermore, WET can be used to charge low power 
devices such as temperature and humidity meters and liquid 
crystal displays 0. Even low-end computation, sensing, and 
communication can be performed by harvesting energy from 
ambient radio frequency (RF) signals including TV, cellular 
networks, and Wi-Fi transmissions 0. 

However, there are several challenges that must be ad¬ 
dressed in order to implement WET Firstly, only a small 
fraction of the energy radiated by an energy transmitter can 
be harvested by the WSN which severely limits the range 
of WET 0, Q. Secondly, the received power levels that 
are suitable for wireless information transfer are not suitable 
for energy transfer, where the absolute received power is of 
interest and not the signal-to-noise ratio (SNR). 
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Massive multiple input multiple output (MIMO) systems, 
where the base station (BS) uses antenna arrays equipped with 
a few hundred antennas, have recently emerged as a leading 
5G wireless communications technology that offer orders of 
magnitude better data rates and energy efficiency than current 
wireless systems 0. Potentially, the use of massive arrays 
could significantly boost the performance of WET as well. 

A. Focus and Contributions 

We consider a scenario where a multi-antenna BS commu¬ 
nicates with and transfers RF power to a WSN. The motivation 
of using an array of antennas is that the BS can exploit 
an array gain, resulting from coherent combination of the 
signals transmitted from each antenna, if it knows the channel 
response. This array gain in turn may increase the operating 
range and/or decrease the amount of transmit energy needed 
to satisfy a given energy harvesting constraint. The drawback 
is that the wireless channel between the BS and the WSN 
fluctuates so that the channel state information (CSI) needs to 
be acquired on a regular basis to enable coherent combining. 

The communication between the array and the sensor takes 
place in two phases as shown in Figure 0 In the first phase, the 
sensor utilizes energy stored in a battery or capacitor to trans¬ 
mit a pilot waveform which is measured at each antenna in the 
BS array, in order to estimate the channel impulse response 
from the sensor. In the second phase, the array beamforms 
energy to the sensor, using the estimated channel responses 
and exploiting reciprocity of the propagation channel. The 
energy harvested by the sensor is used to recharge its capacitor 
or battery, and needed in turn for pilot transmission in phase 
one of the next round and also to perform the main tasks of the 
sensor. In addition, both phases may involve communication 
of information, although that is out of the scope of this paper. 

The main questions asked and answered in this paper are: 
1) What array gain can the massive MIMO setup provide, i.e., 
how does the required uplink pilot energy (and thereby the 
energy storage requirements at the sensor, and the required 
array transmit energy) scale with the number of antennas in 
the array taking into account that all channel responses are 
estimated from pilots? The goal is to determine the scaling 

*We consider time-division duplexing (TDD) mode of communication 
and both the uplink and the downlink communication take place over the 
same narrowband channel. We adopt the widely used reciprocity assumption, 
which implies that the channel gain from the BS to a WSN is the same 
as the channel gain from the WSN to the BS Most physical channels 
satisfy this assumption, but the transceiver hardware might not satisfy this 
condition unless calibration algorithms are applied Go). However, there is 
substantial evidence that such calibration can be performed accurately and 
rather infrequently ED. 
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Fig. 1. Proposed two-phase protocol: Parameters are explained in Section ITT] 

laws for energy transfer using large arrays. 2) How does the 
number of antennas at the BS depend on the path loss or the 
distance between the BS and the WSN? 3) How do the answers 
to the previous questions depend on propagation conditions 
and the correlation between the adjacent antennas? 4) What 
role does power adaptation based on the estimated CSI play 
in improving the outage performance? 

To this end, we derive new expressions that are valid for 
any generic path loss model for the probability of outage 
in energy transfer, defined here as the probability that the 
energy harvested by the WSN is less than the energy that it 
spends on uplink pilots plus the processing energy^ We derive 
expressions for both perfect CSI and imperfect CSI based on 
least squares (LS) or minimum mean square error (MMSE) 
channel estimation and for both Rayleigh fading (without 
dominant channel components) and Rician fading (with domi¬ 
nant channel components). We consider not just the scenarios 
where the downlink array transmit energy is fixed but also 
those where it is adapted based on the channel conditions. We 
present numerical results to quantify the combined effects of 
path loss, energy spent on uplink pilot signaling, the downlink 
energy, the processing energy, the energy harvesting efficiency, 
the Rician if-factor, power adaptation, the correlation between 
adjacent antennas, and imperfect CSI on the probability of 
outage in energy transfer. To summarize, one of the main 
goals of this paper is to estimate the link budget in order to 
determine the feasibility of a system that performs WET using 
multi-antenna arrays. We next discuss the relevant literature on 
WET using multi-antenna arrays. 

B. Related Literature 

The optimal uplink pilot power and the number of antennas 
at the sensor that need to be trained so as to maximize 
the net average harvested energy at the sensor node was 
characterized in lfl2l . However, reference lfl2l did not consider 
the possibility of an outage in energy transfer. The amount 
of time that must be allocated for channel estimation and for 
WET in order to maximize the harvested energy for a multiple 
input single output (MISO) system was investigated in m. 
In m, a wireless powered communication network with one 
multi-antenna BS and a set of single antenna users was studied 
for joint downlink (DL) energy transfer and uplink (UL) infor¬ 
mation transmission via spatial division multiple access. The 
aim was to maximize the minimum data throughput among 
all users by optimizing the DL-UL time allocation, DL energy 

2 We refer to the energy that is needed by the sensor node to perform its 
main tasks as the processing energy. 


beamforming, and UL transmit power allocation. While, the 
optimal training design to maximize the net average harvested 
energy at the sensor over frequency-selective channels was 
studied in 03, energy transfer in the downlink to maximize 
the minimum rate among all users was studied in ES). 

Simultaneous wireless information and power transfer 
(SWIPT), where a multi-antenna BS sends information and 
energy simultaneously to several users which then perform in¬ 
formation decoding or energy harvesting was studied in ca¬ 
ns]]. The authors in lf2Ql investigated when the receiver 
should switch from the information decoding mode to the 
energy harvesting mode based on the instantaneous channel 
and interference conditions so as to achieve various trade-offs 
between wireless information transfer and energy harvesting. 
Receiver design for SWIPT over a point-to-point wireless link 
was investigated in EH- In l22l . the authors studied a hybrid 
network architecture that overlays an uplink cellular network 
with randomly deployed power beacons for charging the mo¬ 
bile devices wirelessly. The tradeoffs between the network pa¬ 
rameters such as transmission powers and the densities of BSs 
and power beacons were derived under an outage constraint 
on the data links. Using a stochastic geometry approach, upper 
bounds on both transmission and power outage probabilities 
for a downlink SWIPT system with ambient RF transmitters 
was developed in |23l . Energy transfer in relay systems to 
simultaneously harvest energy and process information was 
investigated in fl24) - |[26l . 

In contrast to most existing works in the literature, we 
focus on ascertaining whether the use of large antenna arrays 
could substantially extend the feasible range of WET while 
maintaining the receive power level in the same order such 
that a reasonable rectenna efficiency can be maintained. 

The paper is organized as follows: We present the system 
model in Section [Tl] The analysis of the probability of outage 
in energy transfer for different scenarios is given in Section Hill 
and summarized in Tables Q] and [IT] Numerical results and our 
conclusions follow in Section [TV] and Section [VJ respectively. 

The notation X ~ CAM 0. S) means that X is a circularly 
symmetric complex Gaussian RV with zero mean and variance 
<5, and x ~ CAf(m, C) means that x — m is a circularly 
symmetric complex Gaussian random vector with covariance 
matrix C and zero mean vector. The expectation of a RV X 
is denoted by E[2T]. The probability density function (PDF) 
of a RV X is denoted by fx(x). The notation (-)t denotes 
conjugate transpose. Given a complex number z, we denote 
its real part by Re(;z) and imaginary part by Im(;z). 

II. System Model 

We consider a block-fading channel model in which the 
channel impulse response from each antenna at the BS to 
the WSN remains constant during a coherence interval of r 
seconds. The channel realizations are random and they are 
independent across blocks. We, therefore, need to estimate the 
channel after every coherence interval. We assume TDD mode 
of communication so that the channel from the BS to the WSN 
referred to as the downlink channel is the same as the channel 
from the WSN to the BS referred to as the uplink channel. 







3 


Therefore, the BS can take advantage of channel reciprocity 
and make channel measurements using uplink signals. 

We focus on a wireless network where a BS with M 
antennas is used to transfer RF energy to a single antenna 
WSN that has energy harvesting capabilities. We consider a 
scenario where a line-of-sight (LoS) link might be present 
between the BS and the WSN and for which the complex 
channel gain vector h from the BS to the WSN can be 
represented by the Rician fading model as (22 

h = J-^—h d + \! —-— h s , (1) 

where hd £ C Mxl is a deterministic vector containing the line- 
of-sight and the specular components of the channel, f3 denotes 
distance-dependent path loss, K is the Rician factor defined 
as the ratio of the deterministic to the scattered power, and 
h s £ C Mxl denotes the scattered components of the channel 
and is a random vector with i.i.d. zero mean unit variance 
circular symmetric complex Gaussian entries. Furthermore, 
hd = [s/ao s/a ■■■ y/a M -le^w-vW] 1 ' 
where on, i = 0,..., M — 1 denotes the gain of the f h 
antenna which takes a large value if the * th link is good 
and a small value if it is bad, di(4>), i = 1,... ,M — 1 
is the phase shift of the * th antenna with respect to the 
reference antenna and f is the angle of departure/arrival of 
the specular component. Thus, h ~ CM (ji. A/,), where // = 

vW ei( « ••• v^TTe^ ( "- 1)( ^] r 

and Ah = m- By varying K, the model discussed 

above captures a general class of wireless channels spanning 
from a rich-scattering Rayleigh fading channel (K = 0) to a 
completely deterministic channel ( K —> oo). 

A. Uplink Pilot Signaling and Channel Estimation 

The signal y(t) received at the BS when the WSN 
transmits a continuous-time pilot signal y/E u p(t) of duration 
T < t such that J () T |y(£)| 2 dt = 1, is given by 

y(t) = \fEJip(£) + w(t), for t £ [0, T ], (2) 

where E u is the uplink pilot energy in Joule, and h £ C Mxl 
is the channel gain vector from the WSN to the M antennas 
at the BS as defined in <0Q>- Also, w(t) is the thermal noise 
vector at the BS that is independent of h. The objective of the 
pilot signaling is to estimate h given y{t). 

Now, a sufficient statistic for estimating h at the BS is 

y = f p*(t)y(t)dt= + (3) 

Jo 

where w £ C Mxl is the circular symmetric complex additive 
white Gaussian noise (AWGN) at the BS. Furthermore, w ~ 
CM(0, N 0 Im), where N 0 is the noise power spectral density 
in Joule. There are different ways of estimating h depending 
on which type of a priori information that is available at the 
BS. 

3 This is the complex baseband representation of a physical quantity that 
is proportional to the voltage measured across the load connected to the BS 
antenna. The proportionality constant in turn depends on the load resistor 
used. 


1) LS Channel Estimation: This can be used when the 
distributions of the noise and the channel are not known 
a priori. The LS channel estimate is also the maximum 
likelihood estimate in an AWGN setting. Thus, given the 
observation vector y at the BS, the LS channel estimate Jtls 
of h is given by (28l 



This can be simplified to obtain 

h^s = h + JiLS) 


(4) 

(5) 


where Ji LS ~ CM(0, ^-1 m) is the estimation error (28], 
which is circularly symmetric and is independent of h, 
since it is a linear function of w. Furthermore, h^s ~ 

r\r(„ PEu+(.K+l)No f \ 

V yl 1 ’ E U (K+ 1) 1m j ■ 

2) MMSE Channel Estimation: If the distribution of the 
channel and noise are known a priori, MMSE channel estima¬ 
tion can be used. In that case, the MMSE estimate Jimmse of 
h is (28l 


h M MSE = E [h\y] = E [/i]+cov (h, y) (cov (y, y)) 1 (y-E [y]), 

( 6 ) 

where cov (h, y) is the cross-covariance matrix of h and y 
and cov (y, y) is the covariance matrix of y. It is straightfor¬ 
ward to show that cov(h,y) = m and cov(y,y) = 

P E u +(E+i) N o j M ^ Therefore, the MMSE estimate of h in (6]) 
can be simplified to obtain 


ft-MMSE — h + JlMMSE, ( 7 ) 

where Jimmse ~ CAf{0, j 3 e u J 3 (k+i) N o Im ') is the estimation 
error (28l . which is circularly symmetric since it consists of 
noise at the BS and a part of the true channel, both of which 
are circularly symmetric. Also, it is independent of Jimmse- 
Furthermore, h M MSE - CM (/x, (f} EuH K+i)N 0 )(K+ 1 ) 1 ^) • 


B. Transmit Beamforming Based on the Estimated Channel 

In this subsection, we will see how the BS performs transmit 
beamforming based on either the LS or the MMSE channel 
estimate and also characterize the energy harvested. 

1) Transmit Beamforming Based on the LS Channel Es¬ 
timate: Given the channel estimate h^s, the BS performs 
transmit beamforming of energy: it selects the signals emitted 
from the different antennas so that they add up coherently at 
the WSN, i.e., maximizes the harvested received energy. Thus, 

on the downlink, it transmits x(t) = y/Ed ~ LS p' (t), where 

^ LS 

E d is the downlink array transmit energy in Joule, and p (t) 
is a unit energy continuous-time pulse of duration T'. Also, 
T + T' < t. The continuous-time signal y'(t) received by the 
WSN is 

_ , 

y\t) = s/Ed J^ p\t) + w’(t), for t £ [0, T'], (8) 

IIMI 

where w'(t) is the thermal noise at the WSN. Let y denote 
the energy harvesting efficiency of the WSN. Then, the energy 
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harvested Eh in Joule is 


E h 


f]E d 


IIMI 


(9) 


Note that Eh is a random variable since both h and h^s are 
random. We have neglected the contribution from w'(t) to Eh, 
since it is negligible. 

Let us now define 


*LS = 


h[ s h 

fell’ 


( 10 ) 


which is the RV in the harvested energy expression in 
We state below a result that will be used in the performance 
analysis in Section UlI-BI 

Lemma 1: Given the LS channel estimate h^s, 'lb s is a 
complex Gaussian RV with conditional mean 


E 


^LS^LS 


_ PE U _,,r I, 

pE u + {K + l)N 0 ^ Ls|1 

(K + 1)N 0 




and conditional variance 
^LS^LS 


var 


(3E U + (. K + l)JVo ||/i LS ||' 

PN 0 


( 11 ) 


( 12 ) 


pE u + (K + 1)N 0 ' 

Proof: The proof is given in Appendix [A] ■ 

Corollary 1: For Rayleigh fading (K = 0), the RV T IS 
given /ils is distributed as 

PE U ~ pN 0 


H'lsI^ls ~ CM 


PE U + N 0 


IlhLsIl, 


PE U + No 


(13) 


2) Transmit Beamforming Based on the MMSE Channel 
Estimate: If the BS performs transmit beamforming given 
the MMSE estimate /tmmse and on the downlink transmits 
x(t) = sfE) ,~ mm se »'(f) instead, then the signal y'(t) re¬ 


ceived by the 


|| t/,viMsr: 11 

WSN on the downlink is 

t 


s'(i) = + «,'(!). 

11 “'MMSE 11 


(14) 


When using the MMSE estimate for beamforming, the 
energy harvested Eh in Joule is 

2 


E h = r]E d 


^-MMSE^ 


(15) 


||^mmse|| 

We next characterize this RV that is based on the MMSE 
estimate. To that end, let us define 


'I'MMSE — 


h 


MMSE 


h 


(16) 


||^mmse|| 

We state below a result that will be used in the performance 
analysis in Section UlI-CI 

Lemma 2: Given the MMSE channel estimate Jimmse* 
T'mmse is a complex Gaussian RV with conditional mean 


E 


^mmseI^mmse 


= II h 


MMSE | 


(17) 


and conditional variance 

’I'MMSEl/hVIMSE 


var 


PNo 


(18) 


pE u + (K + l)iVo 

Proof: The proof is given in Appendix [B] ■ 

Corollary 2: For Rayleigh fading ( K = 0), the RV Tmmsh 
given /tmmse is distributed as 


’HmmseI^-mmse ~ CM \ ||/»mmse||, ° at 

\ P-t^u H - -/Vq 


(19) 


The conditional statistics derived in this section are used 
subsequently in the analysis of the probability of outage in 
energy transfer in the next section. 


III. Analysis of Probability of Outage in Energy 
Transfer 

Ideally, we want the energy harvested Eh to be greater than 
the sum of the energy E u spent on uplink pilot signaling 
and the processing energy E p that is required by the sensor 
to perform its main tasks. However, this cannot always be 
guaranteed on fading channels. In this section, we compute 
the probability of outage in energy transfer. 

Definition 1: The probability of outage in energy transfer 
P 0 is defined mathematically as 

P Q = Pr (Eh < E u + Ep). (20) 

We compute this probability of outage for scenarios when 
the BS has an LS or an MMSE estimate of the channel 
from itself to the WSN. As a baseline, we also consider the 
case of perfect CSI, in which case the BS knows h exactly. 
This reference case gives us a bound in terms of the best 
outage performance that can be achieved and we include it 
to understand when the uplink pilot is the limiting factor. We 
develop outage expressions not only for scenarios when the 
downlink array transmit energy E d is fixed but also when 
E d is adapted based on the channel conditions. Results for 
different scenarios are summarized in Tables |T] and [TT] 


A. Analysis with Perfect CSI 

As mentioned before, the channel estimation is considered 
error-free if we spend E u on uplink pilot signaling and there 
is no noise in the estimation process. In this subsection, we 
first investigate the scenario where E d is fixed. Thereafter, we 
analyze the probability of outage with power adaptation, where 
Ed is varied based on the instantaneous channel conditions. 

1 ) Without Power Adaptation: With fixed E d , P Q is given 
as follows: 

Theorem 1: For a Rician fading channel, the probability of 
outage I f in energy transfer with perfect CSI and with fixed 
E d is given by 


Po = 1 — Qm 


,\ 


M—1 

2 K J2 a 

2=0 


I2(K + 1)(E U + E P ) 


rjPEd 


where Qm(t) is the ,\/ lh order Marcum-Q function 
Eqn. (4.59)]. 

Proof: The proof is given in Appendix [C] 


(21) 
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TABLE I 

Probability of outage in energy transfer for different scenarios without power adaptation 


Scenario 

Probability of outage in energy transfer 



Perfect CSI, K = 0 

(M-iy. 



LS estimation, K = 0 

1 PE* r -.„( E u +E p \^ M - 1 ( No \ k 

L k ( 

' e u (.e u +e p )\ 

(3E U +N 0 V -nPEd ) 2^k=0 \pE u +N 0 ) 

, vE d N a ) 

MMSE estimation, K = 0 

1 PB " C-' D ( E u +E p \^M- 1 ( No Y 

T . ( 

' E U (E U +E P )\ 

[PE U +N 0 V vPE d ) ^k =0 \f3E u +N 0 J 

P'k \ 

, vE d N a ) 

Perfect CSI, K fO 

i - Qm ( V 2 k E“o 1 °i. V 2 (K+ Z E E : +Ep) ) 



LS estimation, K/0 

No closed-form, can be evaluated numerically using j24) 



MMSE estimation, K A 0 

Single integral form, can be evaluated numerically using |30) 



Note that, if a* = 1, for all i = 0,..., M — 1, then 

Po = 1 ^ Qm {y/2KM, 2{K+1 ^E U d ± ^ 2 ^j ■ Next ’ we state 

the probability of outage in energy transfer for a Rayleigh 
fading channel. 

Corollary 3: For Rayleigh fading ( K = 0) and with fixed 
E d , P 0 is given as follows: 


Po = 1 — Qm 



-y (M 

7 V ’ vP E d ) 

(M- 1)! 


( 22 ) 


where the second equality in (l22l > follows from the identity 
in |29l Eqn (4.71)] and y(-, •) is the lower incomplete Gamma 
function 11301 Eqn. (6.5.2)]. 

2) With Power Adaptation: The probability of outage in 
energy transfer when Ed = | ^ 1 2 , where p is the power control 
parameter, i.e., E d is adapted based on the channel conditions 
is given as follows 0 

Theorem 2: For a Rician or a Rayleigh fading channel, the 
probability of outage P a in energy transfer with perfect CSI 
and with power adaptation can be made zero if and only if 


P > 


E u + E p 
V 


Proof: The proof is given in Appendix [D] 


(23) 


B. Analysis with LS Channel Estimation 

We now investigate the probability of outage in energy 
transfer when the BS performs transmit beamforming using 
the LS channel estimate, first with fixed Ed and thereafter 
with Ed adapted based on the estimated channel conditions. 

1) Without Power Adaptation: With fixed Ed, P 0 for LS 
channel estimation is as follows: 

Theorem 3: For a Rician fading channel, the probability 
of outage P Q in energy transfer with LS channel estimate and 
for a fixed Ed is 


P 0 =£ 


l ~Q 1 vC(^ls),i 


2 (E u +E p fpE u + (K + l)N 0 y 


V/3E d N 0 


(24) 


4 Inverse channel inversion can be impractical for small number of antennas, 
but not as M increases. This is another benefit of having an array of antennas. 


where 


C(^LS) 


2 (^{3E u \\h LS \\ 2 + Re (h^/x) (K + l)iV 0 ) 2 
pN 0 (PE u + (K + l)N 0 )\\h LS \\2 

2Nq(K + l) 2 / Im (^LS^) 

P{PE U + (AT + 1)Nq) I ||hLs|| 

(25) 



Proof: The proof is given in Appendix 10 ■ 

To compute (f24t in closed-form, we need to find the distribu¬ 
tion of C(^ls) given in (f25l) . This is analytically intractable 
but the expectation in (l24l) is easily evaluated numerically. 
A closed-form expression for the outage probability for a 
Rayleigh fading channel can, however, be obtained as stated 
below. 

Corollary 4: For a Rayleigh fading channel (K = 0), the 
probability of outage P Q in energy transfer with LS channel 
estimate and fixed E d is 


P n = 1- 


PE U 


pE u + N, 

M—l 

x E 


k -0 


■ exp 
N n 


PE U + N 0 


E u + E p 

rjpEd 


L k - 


E U (E U + E p ) 
pE d N 0 


(26) 


where hk(-) is the fc th Laguerre polynomial and 


£fc 


1, k < M — 1, 

1 + ^-, k = M — 1. 


(27) 


Proof: The proof is given in Appendix [F] ■ 

Based on Theorems [j] and [3] and Corollaries [3] and |4] 
we observe the following: For fixed M, E u , E p , r/, and 
P, using (l22l > for perfect CSI or using (l26l i for LS channel 
estimation, one can find E d so that a target probability of 
outage in energy transfer is maintained. One can infer how 
the required value of M scales with the path loss B or with 
the distance between the BS and the WSN, for a given P a . 
The loss due to estimation errors can be quantified using the 
analysis in this section. One can also evaluate the role played 
by the LoS component, i.e., the Rician-A" factor on the outage 
probability using (OTb for perfect CSI and using (l24l > for LS 
channel estimation. 
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TABLE II 

Probability of outage in energy transfer for different scenarios with power adaptation 


Scenario 

Probability of outage in energy transfer 

Perfect CSI, K = 0, K f 0 

0 , provided p > 

LS. K = 0 

( 2 V 1 1 x 2 1 l (^h— h— 2 \ ( 1 — x! ^ f /(h—l\ /Zi\\ Xi(^'-xi) 

MMSE estimation, K = 0 

Ef =1 (f )"' 1 e"=o ( 2 T-T 2 ) (-r; 1 ) - 0.(3) - 

LS estimation, K 7 ^ 0 

No closed-form, can be evaluated numerically using ((28) 

MMSE estimation, K f 0 

Single integral form, can be evaluated numerically using d32) 


2) With Power Adaptation: When Ed = -p p is varied 
based on the LS channel estimate, the probability of outage is 
given by the following result. 

Theorem 4: For a Rician fading channel and with power 
adaptation, the probability of outage P a in energy transfer with 
LS channel estimate is 


Theorem 5: For a Rician fading channel (K f 0), the 
probability of outage P Q in energy transfer with MMSE 
channel estimate and fixed Ed is 


Pn = 1- 


2A 1 (K + l) M i 


M -1 


( K zZ o 1 *) 


— exp(-AiA' V 


2=0 


Oi4 


Po = E 


l - Q\ 


C(^-LS), 


2 (E u +E P l(3E u +(K + 1)N 0 ) 


r]/3pN 0 


II^LSI 


(28) 


where C(^ls) is given by (125b . 

Proof: The proof is given in Appendix [G] ■ 

Again, (l28l) cannot be simplified any further but the expec¬ 
tation in |[28l) is easily evaluated numerically. A closed-form 
expression for the outage probability with power adaptation 
and for a Rayleigh fading channel is stated below. 

Corollary 5: For a Rayleigh fading channel and with 
power adaptation, where p > EuEEp ^ , the prob¬ 

ability of outage P 0 in energy transfer with LS channel 
estimation is 


Po = 


M 

E 

ii=i 


2 V 1-1 xl h ~ 2 ^ (2h-l 2 -2\ 
p'j a h-i ) 



where k' = p' = 


ao = 

E u 


2 E u 


/3E U 


N 0 ( f3E u +N 0 )'■ 
ai+bi 


bn = 


2((3E u -\-N 0 ') E u -\-E p 
N 0 r]f3p ’ 


PO = 


pE u +N 0 ’ Ul = 2 a 0 bo 

and X 2 =U 2 ~ y/u\-\. 

Proof: The proof is given in Appendix [H] 


u 2 = 


2po+a 0 -|-b Q 

2aobo 


» Xi = Ul - V U 1 - 


C. Analysis with MMSE Channel Estimation 

In this subsection, we will analyze the probability of outage 
with MMSE channel estimation first with fixed Ed and then 
with Ed adapted based on the estimated channel conditions. 

1 ) Without Power Adaptation: With fixed Ed, P 0 for a 
Rician fading channel is as follows: 


x f y™eM-MK+l)y 2 Q )I M -i 2/Yn K(K + 1) V a it 
Jo \ \ l= o 

X«1 (vS2yo, s j^^p)dyo, (30) 


where Im- i(-, •) is the (M— 1)* order modified Bessel func¬ 
tion of the first kind ||29l Eqn. (4.36)], Q i(-, •) is the first order 
Marcum-Q function |29l Eqn. (4.33)], Ai = /3g "+^ +1)JV ° ; 
and A 2 = 

/V ° n 

Proof: The proof is given in Appendix [T| ■ 

Note that (l30l) is in the form of a single integral in yo and 
probably cannot be simplified any further as the integrand 
involves the product of a modified Bessel function and a 
Marcum-Q function. It is, however, easy to evaluate numeri¬ 
cally. An integral-free closed-form expression for the outage 
probability for a Rayleigh fading channel can be obtained as 
stated below. 

Corollary 6: For a Rayleigh fading channel (K = 0), the 
probability of outage P a in energy transfer for a fixed Ed and 
with MMSE channel estimation is 


P n = 1- 


PE U 


pE u + N, 

M -1 

x E efc 


k -0 


■ exp 
N n 


pE u + N 0 


E u + E p 
ypE d 


L k - 


E U (E U + E p ) 
yE d N 0 


(31) 


where /.;,.(■) is the fc th Laguerre polynomial and Ck is given 
by d27). 

Proof: The proof is given in Appendix Q] ■ 

Note that the expressions for the probability of outage in 
energy transfer is the same for both the MMSE and LS 
estimators under i.i.d. Rayleigh fading and for fixed Ed- 
This is because under i.i.d. Rayleigh fading, the LS and the 
MMSE estimates differ only in terms of the scaling factor. 
Since, we normalize the beamforming vector, this scaling has 
no impact on the end performance. However, under spatial 
correlation or even for uncorrelated Rician fading, MMSE 
gives a better outage performance compared to LS estimator. 
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TABLE III 

Estimates of Path Loss and the Corresponding BS-WSN 
Separation from (7) 


Path loss {ft) 

BS-WSN distance 

60 dB 

7.8 m 

55 dB 

4.1 m 

50 dB 

2.2 m 

45 dB 

1.1 m 


the difference being significant particularly at low SNRs as 
shown in Section lni-D2l and in Figures [8] and [T0]in Section HVl 
2) With Power Adaptation: If, however, Ed = ~ p 

11 h-MMSE 11 

is adapted based on the MMSE estimate, the probability of 
outage in energy transfer for a Rician fading channel is as 
follows: 

Theorem 6: For a Rician fading channel (K f 0) and 
with power adaptation, the probability of outage P D in energy 
transfer with MMSE channel estimate is 


P" = 1 - \m-i - exp(-AiiF <*i) 


j \ > IVl — 1 

( K Ei= 0 a i) 2 


2=0 


x/ 7 ex P( -A 1 (AX-1 ii 1 ,,_ 1 |a\ 


\ 


M—l 

K(K + 1) ^ a iV o 
2=0 


X Qi ^ <E ;„ + Ep) »°) <fr»- (32) 


Proof: The proof is given in Appendix [K] ■ 

Note that (l32l i is in the form of a single integral in yo and 
probably cannot be simplified any further as the integrand 
involves the product of a modified Bessel function and a 
Marcum-Q function. It is, however, easily evaluated numer¬ 
ically. The outage probability for a Rayleigh fading channel 
with power control is given by the following result. 

Corollary 7: For a Rayleigh fading channel (K = 0) and 
with power adaptation where p > Eu + E p , the probability of 
outage P Q in energy transfer with MMSE channel estimate is 


P n = 


M 

E 

n= 1 
2 
2 


n —1 


C 


3n—2 
2 


n— 1 


(! ^ CD 2 - 1 ^ 


i-Ci 

Cl 


n— 1 
c 


2n — c — 2 
n — 1 

Cl(« — Cl) 


i-c? 


(33) 


where k = 


E u +E v 

VP 


p = 


- (.2 

2/3E u [E u +E p __ j 2(pE u +N 0 ) 


N„ 


-,a = 


/3N a 


b = P = «2 = 


_ _ / 2(0E U +N O ) E u +E p n _ 0E u +N o 

V0P _ _ 

2p+ 2 ^ b2 ’ Cl = «i - a/w? - 1, and = v 2 - y/v% - 1. 

Proof: The proof is given in Appendix 0 


Z). Extensions to Non-Identical and Correlated Channels 

1) Independent and Non-Identically Distributed Rayleigh 
Fading Channels: For analytical tractability, we focus on the 
scenario where the M wireless links from the BS to the 
WSN see independent but non-identically distributed (i.n.i.d.) 
Rayleigh fading and derive the outage probability with perfect 
CSI. We do not discuss the case with LS or MMSE channel 


estimation as they are analytically intractable and do not 
provide any additional insights. 

Theorem 7: For i.n.i.d Rayleigh fading channels without 
power adaptation and with perfect CSI, the probability of 
outage P Q in energy transfer is 


Pn = 




M 


M—l 


(i - exp 


3= 1 


Ilk=l,kjtj(^3 Pk) 


(34) 


Proof: The proof involves finding the CDF of |/),|| 2 , 
where each hi ~ CAf (0,/3i) and which can be obtained using 
the result in EH Eqn. (3)]. ■ 

2) Correlated Channels: To analyze the effect of spatial 
correlation on outage probability, we generate the channel 
covariance matrix R using the exponential correlation matrix 
model from lf32l . For this model, the (i,j )th element of R is 
given by 

[®Vki = 



(35) 


This model basically represents a uniform linear array where 
the correlation coefficient between adjacent antennas is given 
by |r|, for 0 < |r| < 1. 

The outage probability analysis developed in this section 
gives us insights about the feasibility of WET using multi¬ 
antenna arrays. 


IV. Numerical Results 

In this section, we present numerical results to quantify 
the potential of using massive antenna arrays for WET using 
the two phase communication scheme in Figure Q] Unless 
mentioned otherwise, we take E u = 10~ 8 J (e.g., 100 pW 
during 100 ps), Ed = 10” 3 J (e.g., 1 W during 1 ms), 
E p = 10~ 7 I (e.g., 1 mW during 100 ps), y = 0.5@ and 

5 The rectenna efficiency depends on the RF input power because the 
rectifier consists of a diode whose impedance varies non-linearly as a function 
of the input power level. The efficiency also depends on the rectenna 
design, the power management unit, the frequency of operation, and the load 
connected to the sensor node a, eu As illustrated in Figure 7 in (4), for a 
specific rectenna design, if the RF power input is above a certain threshold, the 
efficiency does not vary much, whereas below this threshold (which depends 
on the technology used), the efficiency reduces. For the numerical example 
given in Figure 7 of 0, it can be observed that the rectenna efficiency reaches 
almost 45% at an input power level of 1 mW and increases marginally to 50% 
when the input power level is around 6 mW. In our numerical examples, 
we use an efficiency figure of 77 = 0.5 and make sure that the input RF 
power is sufficient to justify this value. For example, with Ed = 10 “ 3 J, 

= —50 dB, and M = 30, the received energy is about 30 x 10 -8 J. 
Over a 100 /i s duration, this corresponds to a received power of 3 mW. So, 
an appropriate justifiable choice of the time duration for the receive signal 
would ensure that the input power is in the order of a mW. As the path loss 
increases or Ed reduces, more antennas can be deployed to ensure that the 
received RF power remains the same. As mentioned above, if the RF power 
received at the sensor is lower than a certain threshold (1 mW for the cited 
rectenna) then the energy harvesting efficiency reduces and a value less than 
0.5 (depending on the input power) will have to be considered (4). Note that 
an intermediate power management unit is needed to impedance match the 
rectenna to the energy storage unit in the sensor node. The design of the 
associated power management unit is discussed in 1341 — [361 . Received power 
levels on the order of a mW are envisioned in applications like charging 
of hearing aids, wireless keyboards and mouses ED As another example, 
the commercially available equipment Powercast (TX91501) harvests about 
4 mW with a 3 W transmitter at a range of 1 meter at 915 MHz earner 
frequency. Another potential application of our work could be in millimeter- 
wave cellular networks where received power on the order of a mW is 
considered reasonable ED. 




































N 0 = k B TlO F / 10 = 10" 20 J, where k B = 1.38 x 10" 23 J/K, 
T = 300 K, and the receiver noise figure is F = 7 dB. 
Also, we consider a uniform linear arrajo for which 0,; (f) = 
2n di cos((/)), i = — 1. We take tf> = 7t/ 3 and 

d = j = 0.06 m, where A is the wavelength at a frequency 
of 2.45 GHz. We also detail below estimation of typical 
product of path loss and energy harvesting efficiency from 
experimental results in 0, that we will use in our numerical 
examples. 

Estimation of path loss and energy harvesting efficiency 
from experiments in IE/.' From a, when a 4 W transmitter 
connected to a vertically polarized fan beam array antenna 
(with gain Gt = 9 dB) is employed, the DC power harvested 
as a function of the distance from the transmitter in a LoS 
situation is plotted in Figure 3 of £7). The carrier frequency 
( f c ) used in the experiment in J7) is 2.45 GHz. Also, this 
experiment was carried out in an office corridor environment. 
From this plot, the product p/3 of the energy harvesting 
efficiency and the path loss can be estimated as follows: 


r)/3 


Pdc 

PtGtGr 


(36) 


where Pdc is the DC power harvested by the sensor, Iffi is 
the transmit power, Gt, and G B are the gains of the transmit 
and receive antennas respectively. We assume that Gr = 0 dB. 
Figure 3 given in 0 shows that at a distance of 1.1 m from the 
transmitter, the DC power harvested is about 0.5 mW. There¬ 
fore, from ( f36b . p/3 = Pt q^q r = 1-58 x 10~ 5 = —48 dB. 
Assuming an energy harvesting efficiency rj = 0.5, this gives 
us a path loss of /3 = —45 dB at a distance of 1.1 m. Similarly, 
different estimates of the path loss and the corresponding 
distances between the BS and the WSN can be obtained as 
listed in Table [III] We vary B around a nominal value of 
—50 dB in our examples. 

Our numerical results are based on analytical expressions 
developed in the paper and since the Monte Carlo simulations 
overlap with these and do not provide any additional informa¬ 
tion, they are not shown @ 

Figure [2] plots P a as a function of M for different values 
of /3 and for K = 2. We observe that by deploying more 
antennas at the BS, a larger path loss (larger distance between 
the BS and the WSN) can be tolerated while keeping the 
outage probability fixed. For example, by going from about 
20 antennas to 100 antennas at the BS, an outage probability 
of 10 6 7 can be maintained even if the path loss increases from 
45 dB to 55 dB. Also, for E u = 10" 8 J and E d = 10" 3 J, 
the difference in performance among MMSE, LS, and perfect 
CSI is negligibly small. 

Figure [3] plots P a as a function of E d for different values 
of M and for K = 2. It can be observed that as E d increases, 
the outage probability decreases. Moreover, as more antennas 
are deployed at the BS, a lower E d is required to keep the 
outage probability at the same value. For example, by going 
from about 10 to 40 antennas at the BS, E d can be reduced 


6 Each antenna in the array is omnidirectional, only the array as a whole 
can form a beam and not each antenna on its own. 

7 Results with LS estimator also overlap with that of the MMSE estimator 

for typical energy levels required to enable WET. 



Fig. 2. MMSE/LS channel estimation: Impact of /3 and M on P 0 (Ed = 
10 “ 3 J, E u = ter 8 J, E p = 10 “ 7 J, No = 10 “ 20 J, K = 2, a t e 1, for 
all i = 0,..., M — 1, and 77 = 0.5). The corresponding perfect CSI results 
are shown using ‘o’. 



Fig. 3. MMSE/LS channel estimation: Impact of Ed and M on P 0 ((3 = 
-50 dB, E u = 10 “ 8 J, E p = 10“ 7 J, N 0 = 1CT 20 J. K 2, a; = 1, for 
all i = 0,.. ., M — 1, and 77 = 0.5). The corresponding perfect CSI results 
are shown using ‘o’. 


by 8 dB, while keeping the outage probability fixed at 10 -6 . 
Thus, the array gain obtained by deploying multiple antennas 
at the BS results in huge savings of radiated energy. One can 
also see that the difference in performance among MMSE, 
LS, and perfect CSI is negligibly small at practical operating 
points. 

Figure [4] plots P 0 as a function of M for three different 
values of the Rician-AT factor, namely, K = 0, K = 2, and 
K = 4 and for both perfect and imperfect CSI obtained again 
using MMSE/LS channel estimation. It can be observed that as 
K increases, the channel becomes more deterministic and the 
outage probability improves with perfect or imperfect CSI. In 
other words, a strong line-of-sight component in the channel 
helps in lowering the outage probability. Also, for the energy 
levels E u = 10” 8 I and E d = HP 3 J, the difference in 
performance among MMSE, LS, and perfect CSI is negligibly 
small. Thus, at these energy levels, one can as well use 
LS channel estimation instead of an MMSE estimator that 
assumes a priori knowledge of the angle of arrival, basically 
h d , the Rician-AT factor, the channel and noise distributions 
without degrading the outage probability relative to the perfect 
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Fig. 4. MMSE/LS channel estimation: Impact of K and M on P 0 ( E u = 
lO " 8 J, Ep = 10 " 7 J, E d = 10 “ 3 I, 3 = -50 dB, N 0 = lO " 20 J, 
a.i = 1, for all i = 0,.. ., M — 1, and rj = 0.5). The corresponding perfect 
CSI results are shown using ‘o’. 



Fig. 5. I.n.i.d channels: Impact of E p and M on P 0 ( E u = 10 —8 J, 
Ed = 10 -3 J, iVo = 10 -20 J, K = 0, Ah = diag(/5i,... ,/3 m), and 
r i = 0.5). Perfect CSI results from oi are shown using different line types, 
corresponding MMSE results using ‘o’ and LS using 

CSI scenario. 

Figure O plots P a as a function of M for different values 
of the processing energy E p and for i.n.i.d. channels. It can 
be observed that by deploying more antennas at the BS, the 
WSN gets higher amount of processing energy to perform its 
main tasks at a given target outage probability. Thus, multiple 
antennas at the BS can help transfer more energy to the energy- 
constrained WSNs. 

As shown in Figure [6j with MMSE/LS channel estimation, 
a too high uplink pilot energy would lead to an increase in the 
outage probability as a larger portion of the harvested energy 
is spent on training. A too low uplink pilot energy would 
again increase the outage probability due to increased channel 
estimation errors. Note that for i.i.d. Rayleigh fading chan¬ 
nels, MMSE and LS estimation give the same performance 
irrespective of E u . With perfect CSI, for smaller E u 's the 
outage performance essentially remains the same, since E u is 
much smaller than E p . 

Figure [7]plots P Q as a function of M for two different values 
of K with and without power adaptation for both LS and 
MMSE estimation. As expected, power adaptation at the BS 



Fig. 6 . Impact of uplink pilot energy on P 0 (E d = 10 3 J, E p = 10 7 J, 
No = 10 -20 J. 0 = -50 dB, K = 0, r =s 0, and ?? = 0.5). 



Fig. 7. Impact of power control on P 0 (No = 10 20 I, 3 = —50 dB, 
Ep = 10 ~ 7 J, E u = 10 -8 J, p = 2.204 X 10 ~ 7 J, m = 1, for all 
% = 0,..., M — 1, and rj = 0.5. MMSE results are shown using linetypes 
and LS results using ‘o’). 

helps to improve the outage probability. For the case when Ed 
is adapted based on the channel conditions, while with perfect 
CSI, the outage probability is zero irrespective of M, with 
MMSE/LS channel estimation it is non-zero but decays very 
quickly to zero as M increases for an appropriately chosen 
value of p for both K = 0 and K = 2. Note that we choose 
p to be greater than or equal to (E u + E p )/p otherwise, the 
outage probability is 1 even with perfect CSI and as stated in 
Theorem 2. 

Figure [8] plots P a as a function of M for perfect CSI, 
LS, and MMSE channel estimation and for two different 
values of the uplink pilot energy E u , basically the low SNR 
regime for uncorrelated Rician fading. Please note that the 
outage performance with perfect CSI is the same for both 
E u = 10 -14 J or E u = 10 -15 J since E u is much smaller 
than E p . While MMSE performs marginally better than LS 
estimation at E u = 10” 14 J, the gap in performance is 
significantly higher at E u = 10 -15 J. 

Figure [9] plots P a as a function of M for different values 
of r. As expected, with an increase in antenna correlation, the 
probability of outage in energy transfer deteriorates. In other 
words, by deploying more antennas at the base station, a higher 
antenna correlation factor can be tolerated while maintaining 
a given target outage probability. Note that for the chosen 
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parameters, MMSE and LS channel estimation perform as well 
as perfect CSI for all values of r. 

Figure [TO] plots P Q as a function of M for perfect CSI, LS, 
and MMSE channel estimation and for two different values of 
the uplink pilot energy E u for spatially correlated Rayleigh 
fading channels. From this figure, it is clear that MMSE 
outperforms LS channel estimation under spatial correlation 
and at low signal-to-noise ratios. While at E u = 10“ 14 J, 
MMSE performs marginally better than LS estimation, the gap 
is significant for E u - 10” 15 J. Note that with perfect CSI, 
the outage performance is the same for both E u = 10” 14 J or 
E u = 10” 15 J, because E u is very small compared to E p . 

V. Conclusions 


Fig. 8. Impact of estimation error on P 0 ( E d = 10 3 J. E p = 10 ' J, 
No = 10 ” 20 J, 0 = -50 dB, K = 2, a ; = 1, for all i = 0,..., M — 1, 
and 77 = 0.5). 



Fig. 9. MMSE/LS channel estimation: Impact of spatial correlation r on P 0 
(E d = 10 ” 3 3, E u = 10 ” 8 J, E p = 10 ” 7 J, N 0 = 10” 2 ° 3,0 = -50 dB. 
K = 0, and 77 = 0.5). The corresponding perfect CSI results are shown using 
‘o’. 


We investigated the feasibility of using multiple antennas at 
the transmitter for WET. Specifically, we derived expressions 
for the outage probability when the BS uses an array of 
antennas to focus and transfer energy to a WSN and where 
the channel from the array to the WSN is estimated using 
pilots sent by the WSN. This is done both with perfect CSI 
and with LS or MMSE channel estimates and for Rayleigh 
fading (without dominant components) and Rician fading 
(with dominant components). We proved that by adding more 
antennas at the BS, we can extend the range for WET while 
maintaining a given target outage probability. We further 
observed that a lower downlink energy is required to get the 
same performance due to huge array gains obtained by multi¬ 
antenna beamforming. 

We observed that for the typical energy levels that are 
used in WET, the difference in outage performance among 
LS, MMSE, and perfect CSI is negligibly small. Further, 
we show that a strong LoS component between the BS and 
the WSN helps improve the outage probability. We also 
show that by deploying more antennas at the BS, a larger 
energy can be transferred to the WSN for it to be able to 
perform its main tasks. While with perfect CSI, outage can 
be completely eliminated by power adaptation based on the 
channel conditions. With power adaptation based on LS or 
MMSE estimate, it can be considerably reduced. 



Appendix 


A. Proof of Lemma [7] 


Conditioned on h^s, T[ s is a complex Gaussian RV with 
mean 


E 


^Lsl^-LS 


U 45 E 

IIMI 


hjh LS 


(37) 


Using standard results on conditional Gaussian RVs (28|, it 
can be shown that 


E 


^I^LS 


M 


I< + 1 




@E U + (A + l)./Vo T ^ 
E U (K + 1) M ) 

x (h LS - p'j . (38) 


Fig. 10. Impact of estimation error under spatial correlation on P Q (Ed = 
10 ” 3 J, Ep = 10 ” 7 J, No = 10” 2 ° 3, 0 = -50 dB, K = 0, r = 0.7 , 
and 77 = 0.5). 
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Substituting E 


in 03 and simplifying yields (ITO . 


Also, the conditional variance is 


var 


'LlsI^ls 


h. 


cov 


(ji\h L s^j h LS 


II^LSII 5 


(39) 


Again, using standard results on conditional Gaus¬ 
sian RVs l28ll . it can be shown that cov = 

p E I m- Therefore, the conditional variance of T'ls 

simplifies to ( 112 b . 


B. Proof of Lemma \2} 

Conditioned on Bmmse, 'I'mmse is a complex Gaussian RV 
with mean 


E 


if 


MMSE | fr-MMSE 


h 


t 

MMSE 


||^mmse|| 


■E 


h\h 


MMSE 


(40) 


Using standard results on conditional Gaussian RVs l l28l . it is 
easy to show that 


E 


h | /tmmse 


= \ B + 

P 2 E U 


P 2 E u 


(pE u + {K+l)N 0 ){K + l) 


L M 


(f3E u + {K + 1 )Nq)(K + 1) 


\ _1 

[m] ( 


^■MMSE — B 


(41) 


Upon simplification, this yields E ^mmseI^mmse 

11 ft-MMSE 11 ■ 

And conditional variance is given by 


’I'mmse I^-mmse 


ft 'MMSE COV ( fc I^MMSE)^MMSE 

j | ^MMSE 11 2 


(42) 


Again using standard results on conditional Gaussian 
RVs ll28l . it can be shown that 

cov(h|h„„ SE ) = — + ^"° + 1 )Wo I„. (43) 

Therefore, the conditional variance of Tmmse simplifies 

to m . 


C. Proof of Theorem [7] 

With perfect CSI, 4' = T'ls = ’Lmmse = = \M, 

where h = h can be either the LS or the MMSE channel 
estimate. Therefore, 


P 0 = PifrjEdl^l 2 <E U + E p ) 
(‘HK + W 

= Pr 


2 (^ + 1 )ii«,i I 2 < 2[K + 1)(E U + E p ) 


P 


V/3E d 


= 1 — Q m 


y 


M—l 


2 K ^2 a i ,' 


l2(K+l)(E u + E p ) 


i—0 


ilPE d 


(44) 


Note that (l44l > follows from the fact that 2 ^ +1 ^ ||/i || 2 is a non¬ 
central chi-square distributed RV with 2 M degrees of freedom 
and non-centrality parameter 2 K 2o 1 a t . 


D. Proof of Theorem [2] 

With power adaptation, E d = Therefore, the proba¬ 

bility of outage in energy transfer is 


Po 


Pr^H/rH 2 


< E u + E p ) — Pr 



E u + E p 
V 



Since the expression is deterministic, P a is one if p < E,,l+Ep . 
If, on the other hand, p > Eu+Rp , then the energy harvested 
Eh = 77 11 /i, 11 2 = rjp > E u + E p and there will never be an 
outage. 


E. Proof of Theorem [7] 

With LS channel estimation, 


Pa =E 


Let T'ls — 
written as 


Pr^LsI 2 < 

1.3 


E u + E v 


h 


LS 


/ Wno 

V 2(4)E„ + (K + 1)V 0 ) 


V E d 

. Therefore, P D in 


(46) 
can be 


Pn = E 


Pr I’J'Lsr < 


2 2 (E u + E P )(/3E U + (.K + 1)N 0 ) 


hhs 


ri(3E d N 0 

(47) 

Using Lemma Q] it can be shown that 
given h,[s, Re('l'is) and lm(T L s) are indepen¬ 
dent Gaussian RVs with conditional statistics 

V2 (0E U 11 h LS 11 2 +Re (hlsIJ.) (K +1)JV 0 ) 

y/f3N o (0E u + (K+l) N 0 ) 11 Sls 11 


E 

E 

var 


Re(4 , Ls)|/t-LS 

Im('h L s)|/rLS 

Re(4' L s)|fiLS 


/2N^(K+l) 


Iml/iJs/i 


= var 


y/P(PE u + (K+\)N 0 ) IlhLsIj 
= 1 . 


and 


Im(4 , L s)|/iLS 

Thus, given h ls, |4 , ls| 2 is a non-central chi-square dis¬ 
tributed RV with 2 degrees of freedom and non-centrality 
parameter given by 


C(^LS) — 


(^f3E u \\h L s\\ 2 + Re (ji LS fij (K + l)No)‘ 
/3N 0 {/3E u + (K+l)N 0 )\\h L s\\ 2 


+ 


2N 0 (K + 1 ) 5 


P{pE u + (I< + 1)N 0 ) 


Im 


V 


II^LsIl 


(48) 


Substituting the cumulative distribution function (CDF) of 
I’LlsI 2 given h^s in (1471) yields (l24l) . 

F. Proof of Corollary 0 

For a Rayleigh fading channel, by substituting K = 0 
in (EU), we get 


P n = E 


1 - Qi 


12 E,, 


PE U 


N 0 pE u + N 0 


II^LSII, 


l2(pE u + N 0 )(E u +E p ) 


V r)/3E d No 

where Qi(-,-) is the first order Marcum-Q function 
Eqn (4.33)]. 


( 49 ) 
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To compute ( l49t . we need to find the distribution of Y = 
II^LsIl = \/I^LSi I 2 4-h I^LSmI 2 - Note that for K = 0, 


hhs 
chi-square 


CM ((), pE £ No ). 1 

ire distributed RV wi 


This implies that 


2 E u y2 


Y 1 is a 


0E u +N o - 

with 2M degrees of freedom since it 
is the sum of the squares of 2 M independent standard normal 
RVs. Therefore, the RV Z = Y 2 has the PDF given by 


To compute d55l ). we need the distribution of Y = ||/i L s|| = 

|/iLSi 1 2 + ''' + I^-lSm 1 2 that we have already evaluated 
in ( |5T1 ). Substituting the PDF of Y from ( |5H in d55l) . we get 


Pn= 1- 


fz(z) = 


Eu. 


M M-l 


exp i Eu \ r ), z> 0 . 


/3E U + N 0 J (M — 1)! ^\pE u + N( 


(50) 


xQi 



» / TP nr 

2M-1 / ^uy 

y exp 


I2E U pE u l2(pE u + N 0 )(E u +E p ) 


N 0 pE u + No 


PP/3N 0 


By transformation of RVs, it can be shown that Y = \[Z = 
11 /tls 11 has the PDF given by 

t ,.s J E u \ M y 2M ~ l ( -E u y 2 \ _ 

/y(^) = 2 , AT 777-TTT ex P oze , . y>°- 



J3E U + N 0 J (M — 1)! ^\I3E U +Nt 


Substituting the PDF of Y from ( 17 1 1 in (l49t . we get 


Note the variable y in both the arguments of the Marcum-Q 
function. If p > Eu ~^ Ep ^ dEu+N 0 \ , <[56b can be simplified 
Eqn. (25)1 tc 


( 51 ) using the identity in ll40l Eqn. (25)] to obtain ( l29l l. 


Pn = 1- 


( Ep > 

\0E u +N o ) 


(M-iy. Jo 


M 

- r>-w ~ Euy2 ) 

Jo y P \pE u + No) 


xQi 


I2E U p E u l2(pE u +No)(E u + E p )\ 


N 0 pE u + N 0 


y , - 


r]/3E d N 0 


Note the variable y in just one of the arguments of the 
Marcum-Q function. Using the identity in (39] Eqn. (9)], 
can be simplified to yield 


dy. 

(52) 


I. Proof of Theorem [5] 

With MMSE channel estimation and for a Rician fading 
channel. 


P n =E 


Pr I^mmseI < 


2 > E u E p 


^MMSE I 


Let 'I'mmse = 


^ MMSE 


/ ft Nq 


pE d 

. Therefore, (l57l) reduces to 


(57) 


2(/3E„+(K + l)N 0 ) 


G. Proof of Theorem 0 

With LS channel estimation and with power adaptation. 


P n = E 


Prfe^< 

Vll^sll 2 ” 


E u + E p 


PP 


Els 


(53) 


Let T'ls — 
written as 

Po 


V 2((JE„ + (E + 1)N 0 ) 


. Therefore, P a in (l53l > can be 


= E 


Pr(JS*!l< 

VH^sll 2 


2(E U +E P XPE U + (K+1)N 0 ) 


p/3pN 0 


h 


’LS 


Po 

= E 


2 2(E u +E p )(l3E u + (K+l)No) 
PrllWMMSE S- 


pf3E d N 0 


RmmseJ 

(58) 


Given h MMSE , Re^MMSs) and Ihi('I'mmse) are independent 
Gaussian RVs. Using Lemma [3 it can be shown that 


E 

E 


Re (’TMMSE ) | ^MMSE 
Im(tE'MMSE) I^MMSE 


\R 


= o. 


PE U +(K+1)N 0 \ U T II 
0Nq ) II"mmse||, 

and the conditional 


variances are given by var 


(54) 


var 


Im (’T MMSE ) | ^MMSE 


Re ( ^MMSE ) | /tMMSE 


From Appendix[E] we know that given h L $, 'Pi_s 1 2 is a non¬ 
central chi-square distributed RV with 2 degrees of freedom 
and non-centrality parameter given by ( l48l ). Substituting the 
CDF of |\Pls 1 2 given fa-Ls in (Edi t yields (l28l >. 


H. Proof of Corollary \5\ 

For a Rayleigh fading channel, by substituting K = 0 
in (l28l> . we get 


= 1. Thus, given Timmse, 

v 1'vimse| 2 is a non-central chi-square distributed RV 
with 2 degrees of freedom and non-centrality parameter 
2 ^ pE 71 +(K+i)No 'j ||^ MMSE || 2 . Therefore, ( l58l ) reduces to 


Po =E 


1-0, (wiiftMMSEii ,( P (i ^7 El,) ' 


where An = 2 ( ^ Eu + (K + 1 l Wo 'l 


(59) 


Po =E 


l — Qi 


12 E u 


P E u 


No pE u + No 


M, 


l2{pE u + N 0 )(E u +E p ) 


pppN 0 


II^LSII 


(55) 


0 _ ^ V p n ° )' 

To compute ( l59l >. we need to find the distribution of 

Yo = ||/immse|| = v/|^mmsei| 2 H-+ |Lmmse m | 2 - R ca n 

be shown that is a non-central chi- 

square distributed RV with 2 M degrees of freedom and non¬ 
centrality parameter 2A 52i=o a t (0E u +(K+i)N o ) . xh ere f 0re . 
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the RV Z 0 = Y 0 2 has the PDF 


Substituting the PDF of Yj from ( f64l > in (l62l >. we get 


fz 0 ( z o) — 


/ 7 -r -t \ M ~l~l 


/r/nr^M-1 \ M. _i. 0 

( K PEi =0 <*i) 2 


M—l 


2=0 


x exp A ( (K + l)^o + K/3 ^ 
x Im-i I 2A, 


xQi 



\ 


M—l 


(65) 


f3K(K + 1) ^ GHzo , > 0, (60) 


2=0 


where A = P e v-+(£^+i)No anc | is the modified Bessel 

function of the (M — l) th order and first kind. 

By transformation of RVs, it can be shown that Yo = 
VZE = II^mmseII has the PDF 


Note the variable yi in just one of the arguments of the 
Marcum-Q function. Using the identity in |39l Eqn. (9)], ( 1651 ) 
can be simplified to yield (OTT) . 

K. Proof of Theorem [7| 

With MMSE channel estimation for a Rician fading channel 
and with power adaptation. 


M—l 


, \ M + 1 

froiyo) = - m~T~ -exp(—AA'/3 a *)2/( 


M—l 

(Kp E 

2=0 


xexp(—A(K+1 )?/q)/m-i 2A 


m 

o 


i=o 


Po =E 

Let 'I'mmse = 


Pr ( l^ MMSE l 2 < Pu + E 1 


\||^mmse|| 2 VP 

-5^“-= and rewrite 


h'MMSE^ 


( 66 ) 


'l ft 5 

V 2(/3E„ + (ir + l)N 0 ) 


\ 


M—l 

@K(K + 1) ^2 a iVo 

i=0 


Po = E 


(61) 


Pr ( I^mmse | 2 < 

\J | h-MMSE 11 2 

2 (E u + E p ) (PE U + (K + l)No) 


Substituting the PDF of Yo from ( 16 1 1 in (l59l > and a simple 
change of variables yields (l30l >. 


J. Proof of Corollary [6] 

For a Rayleigh fading channel, by substituting K = 0 
in (l59b , we get 


vPpNo 


h mmse I 


(67) 

From Appendix |T] we know that given /tmmse, 
( I / mmse| 2 is a non-central chi-square distributed RV 
with 2 degrees of freedom and non-centrality parameter 
2 ^^-HA+iyiVo "j ||fj, MMSE || 2 . Therefore, (l67l) reduces to 


) 


P n =E 


1 - Q\ 


1 2 (PE U + Nq) 


Pn=E 


Wo 


11 ^MMSE 111 


1 — Qi ^V^ollh- MMSeIIj ■ 


I Aq(E u + E p ) | 

VP 


I h 


MMSE | 


( 68 ) 


l2(pE u +N 0 )(E u +E p ) 


r)/3E d N 0 


(62) 


To compute (1621), we need to find the distrib ution of Y\ = 
|I^mmse11 = \/|3immsei| 2 3-+ |^mmse m | 2 - Note that for 


Substituting the PDF of ||/immse|| from (f6Tb in (l68l > and a 
simple change of variables yields l l32l >. 

L. Proof of Corollary [7] 

For a Rayleigh fading channel, by substituting K = 0 
in ( |681 ). we get 


K = 0, h 


MMSE. 


CN (o, ■ 


2 E u \ 


. This implies that p _ g 


0E u +N o J' 

2 C i E'+ N °) yf j s a chi-square distributed RV with 2 M degrees 
of freedom since it is the sum of the squares of 2 M indepen¬ 
dent standard normal RVs. Therefore, the RV Z\ = Yf has 
the PDF (for Z! > 0) 

/ PE U + N 0 \ M zf- 1 f-(PE u +N 0 ) Zl ' 

fz 1 {zi)=[ 09TT , -) — ——exp' 


l-Qi 


JI3E U + N 0 ' UIU 
2 ( PJV 0 111 MMSE|h 


PE U + Nq ^ E u + E p 


Wo 


VP 


| I^MMSE 11 


2 E U J (M — 1)! 


P 2 E U 


(69) 

Substituting the PDF of ||/immse|| from (l64l ) in ( [69b . we get 


(63) 


M 


By transformation of RVs, it can be shown that Y\ = \/~Z~\ = 
| I^MMSE 11 has the PDF (for y ± > 0) 

, , , JPE U + N 0 \ M y™- 1 f — ([3E U + N 0 )y\ 

fyM= 2 {-mr) (M^Tji ex H— we: — 


Po = 1 — 


xQi 


( 0 E u +N 0 \ 

{ 0 2 E U ) 

(M-l)\ J Q 


°°2M-i f-(Wu + N 0 )y 2 
2/i exp 


P 2 E U 


l2(pE u + N 0 )^ I2{(3E U +N 0 )E U + Ep 


Wo 


Wo 


VP 


y i dyi. 


(70) 


(64) 
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Note the variable y \ in both the arguments of the Marcum-Q 

function. For p > Eu+E ‘> , d7 () b can be simplified using the 
identity in J40J Eqn. (25)] to obtain ( |331 ). 
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